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The neutrino self-energy is calculated in a weakly magnetized plasma consists of electrons, protons, 
neutrons and their anti-particles and using this we have calculated the neutrino effective potential 
up to order M^ 4 . In the absence of magnetic field it reduces to the known result. We have 
also calculated explicitly the effective potentials for different backgrounds which may be helpful in 
different environments. By considering the mixing of three active neutrinos in the medium with the 
magnetic field we have derived the survival and conversion probabilities of neutrinos from one flavor 
to another and also the resonance condition is derived. As an application of the above, we considered 
the dense and relativistic plasma of the Gamma-Ray Bursts fireball through which neutrinos of 5-30 
MeV can propagate and depending on the fireball parameters they may oscillate resonantly or non- 
resonantly from one flavor to another. These MeV neutrinos are produced due to stellar collapse or 
merger events which trigger the Gamma-Ray Burst. The fireball itself also produces MeV neutrinos 
due to electron positron annihilation, inverse beta decay and nucleonic bremsstrahlung. Using the 
three neutrino mixing and considering the best fit values of the neutrino parameters, we found that 
electron neutrinos are hard to oscillate to another flavors. On the other hand, the muon neutrinos 
and the tau neutrinos oscillate with equal probability to one another, which depends on the neutrino 
energy, temperature and size of the fireball. Comparison of oscillation probabilities with and without 
magnetic field shows that, they depend on the neutrino energy and also on the size of the fireball. 
By using the resonance condition, we have also estimated the resonance length of the propagating 
neutrinos as well as the baryon content of the fireball. 

I. INTRODUCTION 

The particle propagation in a heat bath with or without magnetic field has attracted much attention due to its 
potential importance in plasma physics, astrophysics and cosmology. The processes which are forbidden in vacuum 
can take place in the medium and even massless particles acquire mass when they propagate through the medium. 
Studying the behavior of particles in such environments requires the technique of thermal field theory. Therefore, in 
connection with these astrophysical and cosmological scenarios it has become increasingly important to understand 
the quantum field theory of elementary processes in the presence of a thermal heat bath. The neutrino self-energy 
is studied in the magnetized medium by many authors, where the effective potential of neutrino is calculated and 
applied in the physics of supernovae, early Universe and physics of Gamma-ray bursts (GRBs)[l|, & d, 0, H, @- 

Gamma-ray bursts (GRBs) are the most luminous objects after the Big Bang in the universe[7[ and believed to 
emit about 10 51 — 10 55 erg in few seconds. During this few seconds non-thermal flashes of about 100 keV to 1-5 MeV 
photons are emitted. The isotropic distribution of GRBsQ in the sky implies that they are of cosmological origin 
0, The GRBs are classified into two categories: short-hard bursts (< 2 s) and long-soft bursts. It is generally 

accepted that long gamma-ray bursts are associated with star forming regions, more specifically related to supernovae 
of type lb and Ic. The observed correlations of the following GRBs with supernovae GRB 980425/SN 1998bw, GRB 
021211/SN 20021t, GRB 030329/SN 2003dh and GRB 0131203/SN 20031w show that long duration GRBs are related 
to the core collapse of massive stars The origin of short-duration bursts are still a mystery, but recently there 
has been tremendous progress due to accurate localization of many short bursts by the Swift [H, [l3| and HETE- 
2[ll| satellites. The afterglow observation of GRB 050709 at z=0.1606[Hj by HETE-2 and the Swift observation 
of afterglow from GRB050709b at z=0.225[]J] and GRB 050724 at z=0.258[lj| seems to support the coalescing of 
compact binaries as the progenitor for the short-hard bursts although definite conclusions can not be drawn at this 
stage. Very rece ntly millisecond magnetars have been considered as possible candidates as the progenitor for the 
short-hard bursts [171 [l8j . For a future study of short-hard GRBs, the ultra-fast flash observatory (UFFO) project is 
proposed[l9j]. 
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Irrespective of the nature of the progenitor or the emission mechanism of the gamma-rays, these huge energies within 
a very small volume imply the formation of and 7 fireball which would expand relativistically. In the standard 
fireball scenario; at the first, a radiation dominated plasma is formed in a compact region with a size cSt ~ 100-1000 
km0, H(j. This creates an opaque 7 — e^ 1 fireball due to the process 7 + 7 — > e + + e~. However, in addition to 7, 
pairs, fireball also contain a small amount of baryons, both from the progenitor and the surrounding medium and 
the electrons associated with the matter (baryons), that increase the opacity and delay the process of emission of 
radiation. The average optical depth of this process is very high. Because of this huge optical depth [2l|, photons can 
not escape freely and even if there are no pairs to begin with, they will form very rapidly and will Compton scatter 
lower energy photons. In the fireball the 7 and pairs will thermalize with a temperature of about 3-10 MeV. 

In this stage, a phase of acceleration begins and the fireball expands relativistically with a large Lorentz factor, 
converting internal energy into bulk kinetic energy. As the fireball shell expands, the baryons will be accelerated by 
radiation pressure. The fireball bulk Lorentz factor increases linearly with radius, until reaching the maximum Lorentz 
factor, so the photon number density and typical energy drop. At certain radius, the photons become optically thin 
(the optical depth is r 77 ~ 1) to both pair production and to Compton scattering off the free electrons associated 
with baryons. At this radius, although much of the initial energy is converted to the kinetic energy of the shell, 
some energy will be radiate away with an approximately black body spectrum. This is the first electromagnetic 
signal detectable from the fireball. For an intermittent central engine with typical variability timescale of St, appears 
adjacent mini-shells with different Lorentz factor, which will collide with each other and will form strong "internal" 
shocks. Later, the fireball shell is eventually decelerated by successive strong external shocks with the ambient medium 
(ISM), propagates into the medium]^]. As in each shell exists a non thermal population of baryons and electrons 
through Fermi acceleration and during each shock the system behaves like an inelastic collision between two or more 
shells converting kinetic energy into internal energy, which is given to the non thermal population of baryons and 
electrons cool via synchrotron emission and/or inverse Compton scattering to produce the observed prompt emission. 
Q. The synchrotron spectrum can be calculate if we know the detailed physical conditions of the radiating region. 
For internal shocks, the so-called 'equipartition' hypothesis is often used, which assumes that the energy is equally 
distributed between protons, electrons and the magnetic field [22|. 

As is well known, during the final stage of the death of a massive star and/or merger of binary stars copious amount 
of neutrinos in the energy range of 5-30 MeV are produced. Some of these objects are possible progenitors of GRBs 
[23j . Apart from the beta decay process there many other processes which are responsible for the production of MeV 
neutrinos in the above scenarios: for example, electron-positron annihilation, nucleonic bremsstrahlung etc, where 
neutrinos of all flavor can be produced 24( . Within the fireball, inverse beta decay as well as electron-positron annihi- 
lation will also produce MeV neutrinos. Many of these neutrinos has been intensively studied in the literature [25II2II] 
and may propagate through the fireball. However, the high-energy neutrinos created by photo-meson production of 
pions in interactions between the fireball gamma-rays and accelerated protons have been studied too[27(. The accre- 
tion disc formed during the collapse or merger is also a potentially important place to produce neutrinos of similar 
energy. Fractions of these neutrinos will propagate through the fireball and they will oscillate [H, [2!| if the accreting 
materials survive for a longer period. Although neutrinos conversions in a polarized medium have been studied [30| . 
resent we have studied the neutrino propagation within the fireball environment with and without magnetic field 
where resonant neutrino oscillate from one flavor to another are studied by considering the mixing of two flavors only. 
In this paper we have calculated the neutrino effective potential in the weak field and done a complete analysis of 
the three neutrino mixing within the magnetized fireball and studied the resonant oscillation of it. By considering 
the best fit neutrino parameters from different experiments, we found that electron neutrino can hardly oscillate to 
other flavor, whereas muon and tau neutrinos can oscillate among themselves with almost equal probability and their 
oscillation probabilities depends on neutrino parameters as well as the fireball parameters. 

The organization of the paper is as follows: In sec. 2, we have derived the neutrino self-energy by using the 
real time formalism of finite temperature field theory [U [32|, HH, H3 and Schwinger's proper-time method (35j |. By 
considering the weak-field approximation we have derived the neutrino effective potential and compare it with the 
effective potential for B = case. We also calculate the effective potential for matter background as well as for 
neutrino background. A brief description about the Gamma-Ray Burst and the fireball model is discussed in sec. 
3. The case of three-neutrino mixing is considered in sec. 4, where we have calculated the survival and conversion 
probabilities of neutrinos and also the resonance condition. In sec. 5, we discuss our results for GRB fireball and a 
brief conclusions is drawn in sec. 6. 



II. NEUTRINO EFFECTIVE POTENTIAL 



As is well known, the particle properties get modified when it is immersed in a heat bath. A massless neutrino 
acquires an effective mass and an effective potential in the medium. The resonant conversion of neutrino from 
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one flavor to another due to the medium effect is important for solar neutrinos which is well known as the MSW 
effect. Similarly the propagation of neutrino in the early universe hot plasma (36j. supernova medium (37j and in the 
GRB fireball [38| can have also many important implications in their respective physics. In all the astrophysical and 
cosmological environment, magnetic field is entangled intrinsically with the matter and it also affect the particle 
properties. Although neutrino can not couple directly to the magnetic field, its effect can be felt through coupling 
to charge particles in the background [39]. Neutrino propagation in a neutron star in the presence of a magnetic field 
and also in the magnetized plasma of the early universe has been studied extensively. 

We use the field theory formalism to study the effect of heat bath and magnetic field on the propagation of 
elementary particles. The effect of magnetic field is taken into account through Schwinger's propertime method[35]. 
The effective potential of a particle is calculated from the real part of its self energy diagram. 

The most general decomposition of the neutrino-self energy in presence of a magnetized medium can be written as: 

E(fc) = R(a0 l{ +a±ft ± + l4 + c^\L, (2.1) 

where fc^ = (k°,k 3 ), kj_ — (fc 1 ,^ 2 ) and u' 1 stands for the 4-velocity of the center-of-mass of the medium given by 
?i M = (1, 0). The projection operators are conventionally defined as R = s(l + 75) and L — |(1 — 75). The effect of 
magnetic field enters through the 4-vector 6 M which is given by 6 M = (0,b). The background classical magnetic field 
vector is along the z-axis and consequently 6 M = (0, 0, 0, 1). So using the four vectors m m and 6 M we can express 

£„ = koi - hi (2.2) 

and the self-energy can be expressed in terms of three independent four- vectors k^_, and 6 M . So we can write 
(£ = RtL) 

£ = a±fc± + bfi + cf (2.3) 

The determinant of jt — S, i.c 

detlfc - t] = 0, (2.4) 

gives the dispersion relation up to leading order in a, b and c as: 

k — |k| = b — c cos <j) — aj_|k| sin 2 (f> — V e ff,B, (2.5) 

for a particle, where <fr is the angle between the neutrino momentum and the magnetic field vector. One has to 
remember that the scalars b and c in this case are not the same if one expresses the self-energy in the form given in 
Eq. (|2.ip . but the V e ff t B is independent of how we express S. Now the Lorentz scalars a, b and c which are functions of 
neutrino energy, momentum and magnetic field can be calculated from the neutrino self-energy due to charge current 
and neutral current interaction of neutrino with the background particles. 

A. Neutrino self-energy 

The one- loop neutrino self-energy in a magnetized medium is comprised of three pieces @, one coming from the 
VF-exchange diagram which we will call £iy(/c), one from the tadpole diagram which we will designate by St(fc) 
and one from the Z-exchange diagram which will be denoted by S^(fc). The total self-energy of the neutrino in a 
magnetized medium then becomes: 

S(fc) = E w (fc) + S z (fc) + E f (fc). (2.6) 



Each of the individual terms appearing in the right-hand side of the above equation can be expressed as in Eq. (|2.1[) 
and the Lorentz scalars a, b and c have contributions from all the three pieces as described above. The individual terms 
on the right hand side of Eq. (12. 6| can be explicitly written as: 



-iE w (k) = J ^(^^^LiS,{p)[^j lv LiW^(q), (2.7) 
zS z (fc)= / £L ( ~ ig \r, L%S v m( k % \ )luLiZ^(q), (2.8) 



(2tt) 4 V\/2cos6W \V2cos9w 



4 



and 



- iZ t (k) = C^-ffl-) Rl, iZ^(0) I y^jTr [ 7 „ (C v + C A75 ) i#(p)] . (2.9) 

The subscripts in E correspond to W-exchange, Z-exchange and Tadpole diagrams. In the above expressions g is the 
weak coupling constant and 6w is the Weinberg angle and g can be expressed in terms of the Fermi coupling constant 
as \[2Gf — g 2 /4:My V - The quantities CV and Ca are the vector and axial- vector coupling constants which come in 
the neutral-current interaction of electrons, protons (p), neutrons (n) and neutrinos with the Z boson. Their forms 
are as follows, 

-± + 2 sin 2 9 W e 
2 — ^ sin (7vf p 



and 



C A = \ "? ^ P • (2.11) 



e, n 



o 



Here H /All/ (g) and stand for the H^-boson propagator and charged lepton propagator respectively in presence of 

a magnetized plasma. The (q) is the Z-boson propagator in vacuum and S vt (p) is the neutrino propagator in a 
thermal bath of neutrinos. The form of the charged lepton propagator in a magnetized medium is given by, 

&(p) = S?(p) + s£(p), (2.12) 

where S®(p) and S^(p) are the charged lepton propagators in presence of an uniform background magnetic field and 
in a magnetized medium respectively. In this article we will always assume that the magnetic field is directed towards 
the z-axis of the coordinate system. With this choice we have, 

/>oo 

iS?(p)= / e^ p ^G(p,s)ds, (2.13) 

where, 

<f>(p,s)=is(pl-rr4- t -^pl). (2.14) 
In the above expression 

P\\ = Pl-Pl, (2-15) 
Pi = P 2 i+Pt, (2.16) 

and z = eBs where e is the magnitude of the electron charge, B is the magnitude of the magnetic held and rag is the 
mass of the charged lepton. In the above equation we have not written another contribution to the phase which is 
e\s\ where e is an infinitesimal quantity. This term renders the s integration convergent. We do not explicitly write 
this term but implicitly we assume the existence of it and it will be written if required. The above equation can also 
be written as: 

<t>(j>, s) = iP{ Po ) - is\p\ + t -^p 2 ± } , (2.17) 

where, 

i>(po) = isipl ~ mj) . (2.18) 

The other term in Eq. (f27T3|) is given as: 

G(p, s) = sec 2 z ifilb + m£(cos 2 z — iT, 3 sin z cos z)] , (2-19) 
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where, 

— Vti - sin 2 z iP ■ u - p ■ b 6 M ) , (2.20) 

= sin z cos z(p ■ u — p ■ b u^) , (2-21) 

and 

S 3 =75^- (2.22) 

The second term on the right-hand side of Eq. (|2.12p denotes the medium contribution to the charged lepton propagator 
and its form is given by: 

/oo 
e^Gip, s) ds , (2.23) 
-oo 

where rjpip ■ u) contains the distribution functions of the particles in the medium and its form is: 

dip ■ u) 9(—p ■ u) 
ripiv ■ u) = —z? — — — r 1 a , , , (2-24 

where [3 and fj,£ are the inverse of the medium temperature and the chemical potential of the charged lepton. 

The form of the VK-propagator in presence of a uniform magnetic field along the z-direction is presented in (39j and 
in this article we only use the linearized (in the magnetic field) form of it. The reason we assume a linearized form of 
the VK-propagator is because the magnitude of the magnetic field we consider is such that eB <C Myj- In this limit 
and in unitary gauge the propagator is given by 



W^iq) 



M 2 W 



3ie 
2M^ 



(2.25) 



where Mw is the W-boson mass. Here we assume that q 2 <C M w and keep terms up to l/M w in the W propagator. 

Let us assume that an electron neutrino v e is propagating in the medium (generalization to other neutrinos is straight 
forward) which contain electrons and positrons, protons, neutrons and all types of neutrinos and anti-neutrinos. 

By evaluating the Eq. (|2.7p explicitly we obtain 



Re S w (fc) = R 



aw±fc ± + bwi- + cw$ 



L, 



(2.26) 



where the Lorentz scalars are given by 



V2G F 



E^iN e - N e ) + k 3 iN° - N°) 



r> poo / 2 TT \ 

■5PJI ^EP-Vo)(--^)(/„. + /.,.) 



(2.27) 



bwo + bw 
V2G F 



1 



3 ml 



2M W 



iN e - N e ) + 



eB , E„ e k 3 



(N° e -N° e ) 



eB 



2tt 2 M w 



/ dp 3 J2V 

J n=0 



6 nfi )i 2k 3 E n 5 nfi + 2E v \E, 



2 E n 



(/e,n ~t~ fe,n) 



(2.28) 



and 



c w 



cwo + cw 



V2Gi 



1 mi 



k 2 



2M W 



M w 



iN»-N») + 



eB 



2n 2 M w Jo 



dp 3Z ^(2-S n>0 ){2E Ue [E n 

n=0 



eB 



KM 



(N e - N e ) 



2E r . 



S n n + 2k 3 E, 



3m| 
2 E„ 



H 

E„ 



ife,n ~t~ fe,n) 



.(2.29) 
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The electron energy in the magnetic field is given by, 



El,n = (Pi +m 2 e + 2neB) = (p§ 



H). 



From Eqs. (f2T28|) and 1(2129]) . we have defined 

coo oo 



6w = -V2d 



eB 



2^M^ 



\ d P3 V (2 - <5„, ) 2k 3 E n 6 n>0 + 2E Ve E n 
J o n=0 I V 



2£L 



(/e,n "i" fe,n) 



and 



1/0 n=0 k v n 



X 5 n ,0 + 2k 3 (E n - ^ - ) }• (/ e , n + fe.n) 



2 E n E n 



In the above equations, the number density of electrons is defined as 

eB 



N e = 



„ 2 X^C 2- <W) / dpzfe^n 

and the number density of electrons for the Lowest Landau (LL) state which corresponds to n = is 

eB f°° 

We can express the Eq. (|2.8[) for Z-exchange as 

ReE z (k) = R(a z ft + b z fi)L, 

and explicit evaluation gives, 

~E Ve ,„ - . 2 1 



V2G/ 



and 



6 Z = V2G f 



{N Ve -N u J 



%Ev 
3M| 



(E„ e )N Ve +(E Ue )N v , 



(2.30) 



(2.31) 



(2.32) 



(2.33) 



(2.34) 



(2.35) 



(2.36) 



(2.37) 



In Eq. (|2.35[) we have a term proportional to ft, because there is no magnetic field. But using the four vectors ji and 
the parallel component of the four vector jt can be decomposed as in Eq. (|2.2p . In the calculation of the potential 

the contribution from these terms will cancel each other and only one which will remain is bz- 
From the tadpole diagram Eq. (|2.9p we get, 



ReE t (k) = V2G F R 



C Vc (N e -N e ) + C Vp (N p - N p ) + C Vn (N n -N n ) + (N Ve - iV„ c ) 



L, 



(2.38) 



So the different contributions to the neutrino self-energy up to order 1/Myy are calculated in a background of e + e , 
nucleons, neutrinos and anti- neutrons. 



B. Weak field limit eB < ml 

In the above subsection, the result obtained is weak compared to the W-boson mass i.e. eB -C M^. But here 
we would like to use another limit eB <C ml that is magnetic field much weaker compared to the one done in the 
above subsection. We also assume that the chemical potential of the background electron gas is much small than 
the electron energy (fi <C E e ). The /i = implies CP symmetric medium where number of electrons equals number 
of positrons. So by taking /i <C E e we assume that N e > N e . In a fireball medium this condition can be satisfied 
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because the excess of electrons will come from the electrons associated with the baryons which will come from the 
central engine. 

In the weak field limit (eB <C ni 2 /e = B c ) and [i <C E e , the electron distribution function can be written as 



(2.39) 



2=0 



Also we shall define 



a = 0fi(l + 1), 



(2.40) 



and 

a = 0m e (l + 1). (2.41) 
Using the above distribution function, the electron number density and other quantities of interest are given below: 



N° - N2 = -4— m 3 



e e 



\ p ..^E(-l)' S inha^(a)] = ^ f |- ] <IV 



* 2 B c 



1=0 



7T 2 \B, 



(2.42) 



N e — N e = ^ V(-l) z sinha 



2 

i=0 



a B c 



(2.43) 



^ jf dp 3 E (f e , + / e , ) = ^ (-) £(-!)' cosha^ (a) 



^(a) 



(2.44) 



2^ / dp 3 — (f e , + / e , ) = ( — ) ^(-1)' cosha tf f», 



'So 



7T 2 \B r ^ 

1=0 



(2.45) 



eB °° f°° m 2 °° 

— -V(2-<5„. ) / dp 3 E n (f e , n + f e . n ) = -f V(-l) z cosha 



(2.46) 



and 



n=0 •'° " ^ 2=0 



-Si (a) 

(7 



eB 
2^2 



■ V(2-5 n , ) / *3^(/e,n + /e,n) = ! ^V(-l) i 

n=0 - 70 ^™ 1=0 



cosh a 



All the above quantities are necessary to evaluate the effective potential. 



I-So(a) 



(2.47) 



(2.48) 



C. Neutrino Potential without Magnetic field 

In the absence of magnetic field the neutrino self-energy and the neutrino effective potential is calculated earlier (4C 
In this case the neutrino self-energy is decomposed as 



ReT,(k) =aft + hi, 



(2.49) 
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and the neutrino effective potential for a massless neutrino is given by 



Veff = b=-±-Tr(Kn< 



(2.50) 



By evaluating the right hand side (RHS) up to order l/Myy gives, 



Also we have 



V eff = V2Gl 



1 



3 m 



2 Ml 



(N e - N e ) B = 



m: 



4 

7T 2 \M W 



E Ve ^(-l)'cosh al 



1=0 



8 \ K x (a) 



(2.51) 



(N e - N e ) B=0 = ?± Y(-\) 1 sinha -K 2 {a). 

7T ^ 



(2.52) 



This is the result obtained in ref.[40| up to order for neutrino propagating in a medium with only electrons 

and positrons in it. 



D. Comparison of V e ff with and without magnetic held 



The neutrino effective potential in a magnetic field is given in Eq. (12. 5p . To simplify our calculation we assume 
that, the magnetic field is along the direction of the neutrino propagation so that cj> = and the aj_ term does not 
contribute. Also one has to remember that by taking B — 0, we should get back the result obtained in Eq. (|2.5ip 
and this is only possible when we take k 3 — ~E U in our calculation. Then the effective potential should be defined as 
(independent of the angle 4> is zero or not), 



V e ff,B = (b~ c)/ k3 =-E„- 

Hence forth we shall replace k% by — E v in our calculation. This gives 



V. 



eff,B 



= V2G F 
eB 

+ 



3 ml 



2 Ml 



eB 

Mr 



(N e - N e ) - 1 + 



eB 



2M 2 , A/ 2 , 



/ d P3 V (2 - S n . Q ) 2E Ue E n S n . Q - 2E Ve 2E n 
J o „=o L V 



771 ; 

E 



(N»-N») 

2 a 

En 



e,n ' J e,n J 



With simplifications this gives, 



(2.53) 



(2.54) 



eff,B 



V2Gi 
B 

2 



1 + 



£(-l)<sinha<| (1 
eB 



1=0 
■2 



3 mi 



2M^ 



m 



2M^ Ml 



,2 \ 2 



E Ve ^(-l)'cosh 



1=0 



ir 2 \M W 

We can write this in a simpler form as 

Vef/,B = 

where the functions $a and $s arc defined as 

oo 

$A = y^(— 1)' sinh a 

1=0 

B 



eB 
M 



5 B 
2~Br 



-K 2 (cr) - J^-K^a] 
a B c 



K (a) +2-4 



B 



V2G F ^ 

TV" 



IV1 W 



l+ 3m 2 



2 Ml 



eB 

Mr 



a 



16 \ K x (a) 
7/i 



2MI 



eB 
Ml 



(2.55) 



(2.56) 



(2.57) 
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FIG. 1: The Eq. (|2.55P is plotted as a function of temperature T/m e for a give B — 01. B c . The unit of V e ff,B is in eV . 
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FIG. 2: Only the magnetic field dependence of the V e //,B is plotted as a function of temperature T/m e for a given B 
The unit of V e ff,B is in eV. 



01.B C 



and 



$s = ^(-l)'cosh 



8 5 B \ . . / B 



16 

7^2 



(2.58) 



By taking B — in Eq. (|2.55p it reduces to Eq. (|2.51|) . So in the weak field limit we get back the potential for B = 
in the medium. Here we have shown only for the W-boson contribution. In Z-exchange diagram we do not have 
magnetic field contribution. In the tadpole diagram only electron loop will be affected by the magnetic field. But as 
the momentum transfer is zero, there will not be any higher order contribution. As the magnetic field is weak, the 
protons and neutrons will not be affected by the magnetic field. 

In Fig. [T] we have plotted the potential Eq. (|2.55|) as a function of temperature in the range to 10 MeV for a fixed 
value of the magnetic field B = 0.1B C . This shows that the potential is an increasing function of temperature. We 
have also shown in Fig. [51 only the magnetic field contribution, i.e. by subtracting the B = part from Eq. (|2.55p . 
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which shows that, the magnetic field contribution is opposite compared to the medium contribution and also order of 
magnitude smaller. 



E. Matter Background 



Let us consider the background with electrons, positrons, protons, neutrons, neutrinos and anti-neutrinos in the 
background. As we are considering the magnetic field to be weak, the magnetic field will have no effect on protons 
and neutrons. For an electron neutrino v e propagating in this background, we have 



V2G F 

"/If 2 

ml v 



E Ve { (N e - N e ) - (7V e ° - N° e ) 
8 \ Ki{a 



(-l)'coshc^ 2- 



B 
K 



(2.59) 



bw + bz + h = b Qe + b 



w 



^ G "[( 1 + ^ + s| + c -) (,v -^- ) + 

+C Vp (N p - N p ) + C Vn (N n - N n ) + 2{N Ve - N Ue ) 

%E Vc 



zB El 



M w M w 



+ (Nu M - N Vll ) + (JV„ T -N Vt )- -jfi [ {E Ve )N Ve + (E V „)N V 



jw, 



(2.60) 



and the coefficient of U is, 



C e = C W + C t — c 0e + c w 

rn 



= V2G t 



eB Ei 



Cw, 



(2.61) 



where bw and cw are given in Eqs. (|2.31[) and (|2.32|) . In the weak field limit these two functions are given as 



>w 



-y/2G f 



2 / m 



2 \ 2 



7T 2 iM, 



^^(-l/cosh 



a 2 2B r 



B r 



a" I a 



(2.62) 



and 



c w 



= V2G F ^ 



2 / mi 



2 \ 2 oo 



7T 2 \M W _ 

Similarly for muon and tau neutrinos, 



E Ve y^(-l)' cosh a 



1=0 



(2.63) 



b M = b Qfl 



V2Gi 



Cy c (JV e - N e ) + C Vp (N p - N p ) + CV„ (JV„ -N n ) + (N Vt - N Ve ) 



+2{N v -N„) + {N Ur -N Ur ) 



8E„ 
3 M 



t((e„Jn^ + (e^)n„ 



(2.64) 



and 



b T = b Qr = y/2Gt 



(N e -N e )+ C Vp (N p - N p ) + C Vn (N n - N n ) + (N Ve - N Ve ) 



+ (N^-N^) + 2(N Ut ~N Vt ) 



8E, 
3M 



2 , {E^)N^ + {E Vt )N^ 



respectively and 



(2.65) 



Cp = co^-CaSK-N^) 
c T = c 0T = -CaSN" - N°). 



(2.66) 
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For muon and tau neutrinos propagating in the medium = c T = 0. The matter potentials experience by different 
neutrinos for cj> = are given as 

V Ve = b e - c e 

V v „ = — Ca c (N® — iV° ) 

V Vr = b T -C Ac (N°-N°). (2.67) 

Putting the values of bi and c; (for I = e, ji, r) one can calculate the neutrino potential in the background. For charge 
neutral matter we should impose 

N e - N e = N p - N p , (2.68) 

and this gives 

C Ve (N e - N e ) + C Vp (N p - N p ) = 0, (2.69) 

in Eqs. (|2.60[> . (|2.64|) and (|2.65p . The particle asymmetry is related to the lepton or baryon asymmetry through the 
relation 

N a ~ N a , 
L a = ^-A (2.70) 



where the number density of photon is 7V 7 = ^C(3)T 3 



F. Only Neutrino Background 

In a newly born neutron star, the neutrinos are trapped because the mean free path of these neutrinos are very 
short compared to the depth of the surrounding medium. So slowly these neutrinos will diffuse out of the region where 
they are trapped called the neutrino sphere. In the neutrino sphere, the different neutrinos have different average 
energy, which are given as fill]: 



(E Ue ) 
(E Vx ) = (E Vm ) 



10 MeV 
15 MeV 
20 MeV, 



(2.71) 



for x = ji,T. If the medium contains only the neutrinos and anti-neutrinos of all flavors, then for propagating v e and 



f M we have 



V2Gf 
8(E, 



3 Mf 



2(N Vc -N Ve ) + (N Vft - #„J + (N Vr - N VT ) 



(E Ve )N VB + (E VB )N Ve 



(2.72) 



and 



y/2Gi 



(N Ve -N Ue ) + 2(N„ -N u ) + (N Ut - N Ur ) 



(2.73) 



respectively and by interchanging \i <-»• t in Eq. (|2.73p we obtain b T for tau neutrino. For only neutrino background 
we have c = 0. As (E Vx ) — (E Va ) in the neutrino sphere and the propagating neutrinos are also in the background, 
in Eq. (|2.73p we take E Vm — (E Ut ) — (E Vj . ) . Now the potential difference between v e and will be 



V efl = b e -b u = V2Gi 



(N Ve -N Ve )-(N V -R v \ 



8 1 



(E„.)[ (E^)N Vc + (E Vt )N Vi 



(E^) 2 N^ + (E u fN v „ 



(2.74) 
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Let us assume that the number density of neutrino and anti-neutrino of all flavors are the same inside the neutrino 
sphere, i.e. 

N Vl =N Vl , l = e,n,T, (2.75) 

and this gives 

V efi = 2.91 x lCT 18 ^, MeV- 2 . (2.76) 

The potential difference between and v T vanishes (V^r = 0) and the potential difference between v e and the sterile 
neutrino v s is given by 

V es =b e -b s = -1.32 x lCT 18 ^, MeV- 2 . (2.77) 

But if we do not take into account the restriction given in Eq. (I2.75P then potential for v ei and v T will different 
from each other. 



III. GRB PHYSICS AND FIREBALL MODEL 



We have already given a short introduction to GRB and fireball model in sec. 1. As we are interested in the 
propagation of neutrinos in the fireball medium, let us discuss a bit more about it. The fireball is formed due 
to the sudden release of huge amount of energy in the form of gamma-rays into a compact region of size cSt and 
it will thermalize with a temperature around 3-10 MeV by producing electron, positron pairs [45|. It will also 
contaminated by baryons both from the progenitor and the surrounding medium which is believed to be in the range 
l(r 8 M Q - 1O" 5 M . 

Among the GRB community, it is strongly believed that the prompt 7 which we see in the rage of few 100 keVs to few 
MeVs for few seconds is due to the synchrotron radiation of charged particles in a magnetic field. But comparatively 
strong magnetic field is needed to fit the observed data. But it is difficult to estimate the strength of the magnetic 
field from the first principle. One would expect large magnetic field if the progenitors are highly magnetized, for 
example, magnetars with B ~ 10 16 G. A relatively small pre-existing magnetic field can be amplified due to turbulent 
dynamo mechanism, compression or shearing. Also under suitable condition the neutrino-electron interaction in the 
fireball plasma will be able to amplify pre-existing small scale magnetic field. Despite all these, there is no satisfactory 
explanation for the existence of strong magnetic field in the fireball. Also even if some magnetic flux is carried by 
the outflow, it will decrease due to the expansion of the fireball at a larger radius. But the strength of the magnetic 
field will be much smaller than the critical field Be- So the derivation of the effective potential for weak field limit is 
justified here. However, if we can measure the polarization of the GRBs, it will be helpful to estimate the magnetic 
field in the fireball as well as give information about the nature of the central engine. 

Here we consider a CP-asymmetric 7 and e~e + fireball, where the excess of electrons come from the electrons 
associated with the baryons within the fireball. Here for simplicity we assume that the fireball is charge neutral 
L e = L p and spherical with an initial radius R ~ (100 — 1000) km and it has equal number of protons and neutrons. 
Then the baryon load in the fireball can be given by 

M b ~ ^( 3 )L e T 3 R 3 m p 

~ 2.23 x IQ-^LeTl^yR^MQ. (3.1) 

where Tm & v is the fireball temperature expressed in MeV and lies in the range 3-10. The quantity Rj is in units of 
10 7 cm and m p is the proton mass. For ultra relativistic expansion of the fireball, we assume the baryon load in it 
to be in the range 10 -8 Mq — 1O~ 5 M0 which corresponds to lepton asymmetry in the range 8.1 x 10~ 4 i?f 3 < L e < 
8.1 x lO^flf 3 . 

We have already discussed about the different origins of 5-30 MeV neutrinos. Once these neutrinos are produced, 
fractions of these neutrinos may propagate through the fireball which is in an extreme condition and may affect the 
propagation of these neutrinos through it. 



IV. THREE-NEUTRINO MIXING 

To find the neutrino oscillation probabilities, we have to solve the Schroedinger's equation, given by 

(4.1) 



13 



and the state vector in the flavor basis is defined as 

v={v e ,v^v T ) T . (4.2) 

The effective Hamiltonian is 

H = U ■ H$ -W + diag(V e , 0,0), (4.3) 

with 

ffo = ^diag(-Am 2 21 , 0, Aj 2 ). (4.4) 
Here V e is the charge current (CC) matter potential and U is the three neutrino mixing matrix given bv[43l [44| 

:) 



/ C13C12 S12C13 Si 3 \ 

U = I -S12C23 - S23S13C12 C23C12 - S23S13S12 S23C13 ] , (4.5) 

V S23S12 - S 13 C 2 3Ci2 -S23C12 - S13S12C23 C23C13 , 



where sy = sin 6V, and Cy = cos Oij. For anti-neutrinos one has to replace (N a — N a ) by — (N a — N a ) and U by U*. 
The higher order contribution to potential does not change the sign. Here we have to emphasize that the neutral 
current (NC) contribution is not taken into account. This is because in the matter background the NC contribution 
to all the neutrinos is the same and when we take the difference of potential, this contribution will be cancelled out 
and does not affect the neutrino oscillation. But it has to be remembered that, in the neutrino background where 
N v — 7^ 0, the potential for different neutrinos are different which described in Sec lII Fl and in this case we can 
not neglect the NC contribution. The different neutrino probabilities are given as 

Pee = 1 — ^ s 13,m C 13,m^31' 

P HH = 1 ~ 4s 13,m c 13,ro s 23'S , 31 ~ 4 «13,m s 23 c 23'$81 ~ 4c 13,m s 23 c 23^32, 

D 1 A 2 2 4o A 2 22o A 2 22c 

^tt - 1 - 4s 13,m c 13,m c 23*31 - 4S 13 m S 2 3 c 23^21 - 4C 13 m S 2 3C 23 032, 

Per = 4s 13 m C 13 m C 23 S'31 



(4.6) 



where 



and 



sin2fl 13 m = Sin2gl3 (4.7) 

V(cos20 13 - 2£„V r e /Am§ 2 )2 + (sin20 13 ) 2 



2 \sm20i3 )5 

2 ^ m 32 / sm 2^13 



where 



AM32 = — ^ • 9fl +1 +^Ve 

« 9 a 9 sin 26*13 / , „s 

A/4 = Am 2 32 , 4.9 
sm26»i 3im 



sin 2 6>i 3 ,m = 2 I 1 \/"' -*m-Mt-:.. 



1 



cos' 13 , m = - ( 1 + \J 1 - sin 2 20 13>m ) (4.10) 
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TABLE I: We have shown the different observables of the fireball as well as the neutrino resonance length L res for Am| 3 = 10~ 2S, eV 2 
and B = O.LB c . 



E v ,Me.V 


T(MeV) 


li(eV) 


L e 


Lres(cm) 


M b (R$M Q ) 


5 


3 


2.47412 


1.099 x lO - " 


4.737 x 10 b 


7.064 x 10 _a 




10 


0.960329 


1.284 x 10 -7 




3.057 x 10" 8 


10 


3 


1.33789 


5.943 x 10"' 


9.474 x 10 b 


3.819 x 10 _y 




10 


1.59669 


2.135 x 10" 7 




5.083 x 10" 8 


20 


3 


0.869801 


3.864 x 10 _Y 


1.895 x 10' 


2.483 x 10" M 




10 


3.03194 


4.055 x 10" 7 




9.653 x 10" 8 


30 


3 


0.804489 


3.574 x 10 _Y 


2.842 x 10 Y 


2.297 x 10" y 




10 


4.50277 


6.022 x 10" 7 




1.434 x 10" 7 



TABLE II: We have shown the different observables of the fireball as well as the neutrino resonance length L res for Am.23 = 

W- 2 - 2 eV 2 and B = 0.1B C . 





T(MeV) 


fi(eV) 


L e 


L res (cm) 


M b (R$M @ ) 


5 


3 


12.1177 


5.383 x 10" 


a 


9.452 x W b 


3.460 x 10" s 




10 


1.82987 


2.447 x 10" 


7 




5.826 x 10" 8 


10 


3 


6.17192 


2.742 x 10" 


6 


1.890 x 10° 


1.762 x 10" 8 




10 


2.03432 


2.721 x 10" 


7 




6.477 x 10 -8 


20 


3 


3.27657 


1.456 x 10" 


6 


3.780 x 10 B 


9.355 x 10 _y 




10 


3.24878 


4.346 x 10" 


7 




1.034 x 10" 7 


30 


3 


2.4178 


1.074 x 10" 


6 


5.671 x 10° 


6.903 x 10 _y 




10 


4.6475 


6.216 x 10" 


7 




1.480 x 10" 7 



The oscillation length for the neutrino is given by 



L osc = . = (4.11) 

cos 2 20i3(l - a l E - v ' ) 2 + sin 2 26»i3 



Am^ 2 cos 2^13 

where L v = AttE v / 'Am§ 2 is the vacuum oscillation length. For resonance to occur, we should have V e jf,B — V e > 
and 

2E V V K 



J v v e 
'32 

By putting V e and simplifying we obtain 



cos2# 13 = — A^. (4.12) 
Ami 



$ A - 1.58027 X 10- w E u>M eV$B ^ 2.24208^^_cos26»i3, (4.13) 

where Am§ 2 is expressed in eV 2 and E u _MeV is in MeV. The functions <&a and <£>b are defined in Eqs. (|2.57[) and 
(|2.58p . Now we have to evaluate the above condition for given values of Ato 2 2 and cos26*i3 from experiments and 
different values of temperature (T) and chemical potential (fi). At resonance, the oscillation length becomes the 
resonance length and can be given by 

(4.14) 



sin 26 



13 



So far we have assumed that the neutrino potential does not vary with distance. However V e ff t B will vary with 
distance. So we have to consider the adiabatic condition at the resonance, which can be given by 

2 ^ m K^\U d ^Y >i 



7T V 2E V J V dr 



2 



2 



3.62 x 1(T 2 ( -^2_ sin 2 6 13 ) — > 1, (4.15) 
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FIG. 3: The contour plot of the resonance condition of Eq. (|4. 13(1 as functions of T/m e and log [^-1 is shown for different neutrino 
energies and B = 0.1B C where (a) is for Am§ 2 = 10" 2 ' 9 eV 2 and (b) is for Am| 2 = 10" 2 ' 2 eV 2 . 
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FIG. 4: The survival probability of muon neutrinos P^p, is plotted as a function of log 
(a) and L = 1000 km (b). The neutrino energy and magnetic field are shown in it. 



eV 2 



for the fireball radius L = 100 km 



where 



$' = ^- 1.58027 xlO-^w^ 



dx : dx 

In the above equations we have expressed l cm in centimeter and a; is a dimensionless variable. 



(4.16) 
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log[A n»2/eV A 2] 



FIG. 5: The probability P efl is plotted as a function of log 
The neutrino energy and magnetic field are shown in it. 



-2.9 -2.8 -2.7 -2.6 -2.5 -2.4 -2.3 -2.2 
log[Anf2/eV A 2] 

, for the fireball radius L = 100 km (a) and L = 1000 km (b) 
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FIG. 6: This is same as Fig. 5 but for neutrino energy E v = 30 MeV 



RESULTS 



We have done a complete analysis for three neutrino mixing with and without magnetic field. For our analysis we 
have used the result given in ref.[43|, 



1.4 x 1CT 3 < Amj 2 /eV 2 < 6.0 x 10" 
#13 s 6° 

32° < 6> 23 < 60°. 



(5.1) 



The above result is obtained by performing a global analysis and taking full set of data from solar, atmospheric and 
reactor experiments. In the above we consider 623 — 45°. 

Different values of /i and T are shown in Fig. 3 for which the resonance condition in Eq. (14.13|) is satisfied. We 
have used two extreme values of Am| 2 i-e. 10 -2 9 eV 2 in Fig. 3a and 10~ 2 21 eV 2 in Fig. 3b for B = 0.1B C , #13 = 6° 
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and four different neutrino energies 5, 10, 20 and 30 MeV respectively. In both small and large values of Amfj 2 and 
fireball temperature in the range 3 to 10 MeV, the chemical potential of the electron is in the 1-12 eV range which 
are shown in Table I and II. In Table I we have shown the different fireball observables for Am§ 2 = lO~ 29 eV 2 . It 
shows that going from neutrino energy 5 MeV to 30 MeV, the resonance length changes between 47 km to 284 km 
and lepton asymmetry ~ 10 -7 to 10~ 6 . For a charge neutral plasma L e = L p this translates to a baryon load in 
the fireball in the range 2.3 x 10 _9 M Q < < 1.4 x 10 _7 M Q . In Table II we have done the same analysis but for 
Am§ 2 = 10~ 22 eV 2 . Here the chemical potential is higher compared to the one in Table I. This shows that shift in 
Am| 2 towards higher value also shift the \i in the same direction. In this case there is not much change in L e and Mj, 
and the resonance length lies in the range 9.4 km < L res < 57 km. 

The above analysis of the resonance condition shows that, in the temperature range of 3 to 10 MeV, the resonance 
condition is satisfied for electron chemical potential (fi) between 1 to 12 eV. Also for neutrino energy in the range 5 
to 30 MeV, the resonance length lies below 284 km, which shows that neutrinos can resonantly oscillate within the 
fireball of radius 100 to 1000 km. The baryon load of the fireball also lies in the range 10" 9 M© < M b < 10~ 7 M Q . If 
the potential has a profile that means, the functions &a and $b depend on the length scale then we found that for 
10 -2.9 e y2 < Am | 2 < iq-2.2 e y2 and 5 < E w . Me v < 30, the parameter l cm /& will lie in the range 10~ 10 to 2.5 x 10~ 9 
to satisfy the condition given in Eq. (|4.15[) . 

The survival and conversion probabilities for the active neutrinos are plotted as function of Am 2 2 in the range 
IQ-2.9 eV 2 tQ 10 -2.2 eV 2 for B = and B = q IBc for twQ neutrino energies 5 MeV and 30 MeV from Fig. 4 to 

Fig. 7. We have also considered two different length scales for the fireball i.e. 100 km and 1000 km to see how the 
probabilities changes when the length scale of the fireball changes. As we are taking #23 = 45°, the probabilities 
P ef i = Per and also P w = P TT . 

We have plotted the survival probability of muon neutrino for L — 100 km in Fig. 4a and for L — 1000 km in Fig. 
4b. The survival probability of muon P MM neutrino in Fig. 3a, for neutrino energy 5 MeV and magnetic field B = 
is 180° out of phase compared to the same neutrino energy but for B = 0.1B C . For B = case the probability varies 
between 0.6 and unity and for B = 0.1B C it is between and 0.4. Going from 5 MeV to 30 MeV we saw that, for 
B = 0,the P^fj, decreases and lies between 0.3 and 0.5 and for B = Q.1B C lies between 0.42 and 0.68. Going from 
L = 100 km to L = 1000 km (Fig. 4b), we see that both B = and B = 0.1B C have a small phase difference and the 
Pfj,fj, varies between and 1. But P MM for B = 0.1£> c lags behind the one for B = for both neutrino energies 5 MeV 
and 30 MeV. The probability for low energy (5 MeV curve) neutrino oscillates faster than the one for 30 MeV. 

In Fig. 5, we have plotted the conversion probability P efJl for E u — 5 MeV, L = 100 km (5a) and E v = 5 MeV, 
L = 1000 km (5b) respectively for both B — and B = 0.1B C . It shows that in Fig. 5a, the P efl for B = and 
B = Q.1B C are having the same phase difference of 180° as in Fig. 4a. But the probability is very small ~ 10~ 9 . 
Going from Fig. 5a to 5b (L — 1000 km) , we saw that the phase difference is almost gone away and P efl oscillates 
much faster compared to the one in Fig 5a. 
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In Figs. 6a and 6b, we have the same probability P efl as in Fig. 5, but here the neutrino energy E v = 30AfeV '. It 
is clearly seen in Fig. 6a that both the probabilities are out of phase and are very small ~ 1CP 10 . On the other hand 
for L = 1000 km (Fig. 6b) the phase difference is gone and the P efJ , oscillates much faster than the one in Fig. 6a. 
We have the P efl = P eT ~ 10~ 10 which gives P ee ~ 1. This shows that the electron neutrinos propagating within the 
fireball can not oscillate to other neutrinos. 

In Fig. 7a and 7b we have plotted the P^ T for L = 100 km and L = 1000 km respectively. For L = 100 km the 
probability for B = and B = 0.1-B C are out of phase for both E y = 5 and 30MeV in Fig 7a. In Fig. 7b for 
L = 1000 km there is a small phase difference between the B = and B = 0.1i? c probabilities. Comparison of Fig. 
4a with Fig. 7a and Fig. 4b with Fig. 7b show that the B = probability (in Fig. 4a and Fig. 7a) and B — 0.1B C 
probability (in Fig. 4b and Fig. 7b) are 180° out of phase. We obtain this because the probability satisfies the 
condition 

P^+Pe„+P„ T = 1- (5.2) 

We have shown in Figs. 5 and 6 that P eM = P eT and they are very small which gives P^ T ~ 1 — P w . The P TT is same 
as Pfifi - 

From our analysis we see that P ee ~ 1 and is almost independent of the energy of the neutrinos and the size of 
the fireball, which shows that for small mixing corresponding to 6*13 = 6° and for 823 = 45°, the electron neutrino 
almost does not oscillate to any other flavor which is obvious from the Fig. [5] and Fig. [6j On the other hand, the 
muon and tau neutrinos oscillate among themselves with equal probability and the oscillation depends on the neutrino 
energy, magnetic field and size of the fireball. Comparison of B = and B 7^ results show that the magnetic field 
contribution is order of magnitude smaller than the medium case. But depending on the size of the fireball, the 
probability for B — and are either in phase or out of phase. 

VI. CONCLUSIONS 

We have shown that neutrino self-energy in the presence of a magnetic field can also be expressed as 

£ = a ± ft ± +bii + df>, (6.1) 

by absorbing the fit, component with the two four vectors jh and p. The above decomposition is only valid when 
the magnetic field is along the z-axis. In addition to this we have also shown that the neutrino effective potential 
Veff,B is independent of how we decompose the the self-energy in terms of Lorentz scalars as shown in Eqs. (|2.1j) and 
(|6.ip . We have explicitly calculated the Ve//,B up to order M^ 4 in the weak field limit eB <C ml in terms of Bessel 
Functions and recover the result for B = limit which can only be obtained when £3 the third component of the 
neutrino momentum is replaced by —E v . We have also calculated the neutrino effective potential when the background 
contains (i) e~,e + , protons, neutrons and neutrinos and (ii) only neutrinos in the background. By considering the 
three-neutrino mixing we have studied the active-active neutrino oscillation process v a <-* Uf, (a and b are active) in 
the weakly magnetized e~e + , p and n plasma of the GRB fireball assuming it to be spherical with a radius of 100 to 
1000 km and temperature in the range 3-10 MeV. We further assume that the fireball is charge neutral due to the 
presence of protons and their accompanying electrons. The baryon load of the fireball is solely due to the presence of 
almost equal number of protons and neutrons in it. 

Our analysis shows that the v e almost does not oscillate to any other flavors and P ee — 1 is independent of the v e 
energy as well as the background magnetic field. The non oscillation of v e to other flavors gives P efi and P eT very small 
of the order of 10~ 10 . But the and v T oscillate among themselves, which depends on the energy of the neutrino, 
magnetic field and also on the size of the fireball. We analyzed our result by taking two different radius of the fireball 
i.e. 100 km and 1000 km, neutrino energy in the range 5 to 30 MeV and magnetic field B = and B — Q.1B C with 
B c = 4.14 x 10 13 Gauss. We found that the probability for B = and B = 0.1B C are out of phase by 180° for L = 100 
km and almost in phase for L — 1000 km. For L = 100 km, the P MM and P^r do vary between 0-0.5 or 0.5-1. On the 
other hand for L = 1000 km, it varies between and 1. Also in this case low energy neutrinos oscillate faster than 
the high energy one. 

We have also analyzed the resonance condition and found that, to satisfy the resonance condition, the electron 
chemical potential in the fireball lies in the range 1-12 eV. For neutrino energy in the range 5 to 30 MeV the 
resonance length lies in the range 9.4 km to 284 km. So if we consider a fireball of 100 km to 1000 km radius this 
shows that the resonant oscillation of and v T neutrinos can take place but not the v e . The baryon load calculated 
by using the resonance condition lies in the range 1O _9 M0 < Mf, < 10~ 7 M©. Depending on the size of the fireball 
the probability for B = and B ^ are either in phase or out of phase. 
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